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Abstract 

This work investigates the collective motion of kaons in heavy ion reactions at SIS energies 
(about 1-2 GeV/nucleon). A radial collective flow of K + mesons is predicted to exist in 
central Au + Au collisions, which manifests in a characteristic "shoulder-arm" shape of 
the transverse mass spectrum of the midrapidity K + mesons. The K + radial flow arises 
from the repulsive K + mean field in nuclear matter. In spite of a strong reabsorption and 
rescattering the attractive K~ mean field leads as well to a collective radial flow of K~ 
mesons. The K~ radial flow, however, is different from that of K + mesons and can be 
observed by a characteristic " concave" structure of the transverse mass spectrum of the K ~ 
mesons emitted at midrapidity. The kaon radial flows can therefore serve as a novel tool for 
the investigation of kaon properties in dense nuclear matter. 
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I. INTRODUCTION 



Collective motion in heavy ion reactions has attracted increasing experimental and theo- 
retical interests. "Bounce-off" |T]-^|, "squeeze-out" @,||] as well as radial flow of nuclear 
matter is predicted by theories and observed in experiments. By "bounce-off" one means 
a deflection of the spectators sidewards in the reaction plane, while "squeeze-out" denotes 
nuclear matter jets in the direction perpendicular to the reaction plane emitted from the 
high density overlap region ( so-called fireball). Radial flow is a collective expansion of 
the compressed nuclear matter. The collective motion has proved to be a good probe for 
the equation of state of nuclear matter (EOS) fl9j|i"0| as well as for the heavy ion reaction 
dynamics ||l"l~Hl3ll- 



The techniques developed in studying the collective flow of nuclear matter have also been 



applied to secondary particles such as pions, kaons, hyperons etc. [[L4j -|26|] . For secondary 
particles "squeeze-out" and "transverse flow" are defined in the same way as for nucleons 
by the azimuthal anisotropy and the rapidity dependence of the in-plane component of the 
particle momentum, respectively. The prediction of a pion transverse flow anticorrelated 
with that of nucleons has been confirmed by experiments [|14] --|T6|| . Kaon transverse flow 



attracts also attention since it might provide information about possible modifications of 
kaon properties in dense nuclear matter flT7|-|T9[|. The in-medium kaon properties are an 
important issue not only for nuclear physics, but also for astrophysics with respect to the 
evolution of compact stars E7[. Secondary particles are usually created from elementary 



hadron-hadron collisions which take place randomly in the colliding system. In the spirit of a 



proposal by Fermi ||28|| , multiple particle productions from energetic hadron-hadron collisions 
can be understood in terms of a thermal description. Due to the fairly strong interaction 
between particles involved in a collision, the available phase space will be occupied according 
to statistical laws. This physical picture has been supported by extensive experiments on 
particle production from heavy ion reactions. Therefore, a collective motion of secondary 
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particles should be identified by observing deviations from the pure thermal description. It 
still remains an open question if secondary particles can exhibit collective motion in heavy 
ion reactions in this strict sense. 

In this paper, we explore the possible existence of a collective motion of K + and K~ 
mesons in heavy ion reactions by comparing to a pure thermal picture. We mainly consider 
Au + Au reactions at SIS energies ( about 1-2 GeV/nucleon ), since this system has been 
extensively studied by the experimentalists at GSI/Darmstadt. Secondly, at an incident en- 
ergy below the the kaon production threshold in free NN collisions ( about 1.58 GeV/nucleon 
for K + mesons and 2.5 GeV/nucleon for K~ mesons), one can expect a low thermal velocity 
of the kaons. Consequently, a collective motion, if it exists, will show up clearly. In the 
present study we are especially interested in kaons emitted at midrapidity. These kaons are 
most likely produced in compressed nuclear matter stopped in the overlap region of colliding 
nuclei ( the so-called fireball ). The reaction dynamics is thereby described within the trans- 
port approach of Quantum Molecular Dynamics (QMD). We find in this work that both the 
midrapidity K + and K~ mesons exhibit a novel collective motion in radial direction, which 
we call kaon radial flow. We demonstrate that the kaon radial flow can provide information 
on in-medium modifications of kaon properties. 

This paper is organized as follows. In Section 2, we briefly sketch the QMD approach 
and describe the treatment of the kaon dynamics. In Section 3 we elaborate the signature 
of the K + radial flow and demonstrate its collective nature. In Section 4, we study the K~ 
radial flow, which is found to be quite different from the K + flow. The role of rescattering 
and reabsorption of the K~ mesons is also investigated. In Section 5 we summarize the 
main results and conclusions of this work. 

II. DESCRIPTION OF KAONS IN QMD 
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A. QMD model 



The QMD model [29] has been extensively used for the investigation of heavy ion re- 
actions. In addition to nucleons, we include also A and N* resonances, which turned out 
to play a non trivial role at an incident energy of 1-2 GeV/nucleon PU||?T|| . Pions are pro- 
duced by the decay of these resonances. Their Coulomb interaction with the baryons is also 
included. This is important for a correct description of the pion dynamics, since the pion 
phase space is highly isospin dependence [B2,B3,B3[. 



B. Sources of K + and K mesons 



In heavy ion reactions, K + and K~ mesons are created from elementary hadron-hadron 
collisions. Since only very energetic hadron-hadron collisions are able to produce a kaon 
(or antikaon), one usually assumes that no coherence exists between the different collisions 
producing kaons. At SIS energies the major processes responsible for kaon and antikaon 
production are baryon-baryon collisions 

B + B -> B + Y + K + (1) 

B + B ^ B + B + K + (K°) + K- (2) 

and pion-baryon collisions 

n + B^Y + K + (3) 



tt + B B + K + (K°) + R-. (4) 

Here B stands for a nucleon or nucleonic resonance while Y denotes a hyperon ( £ or A ). At 
SIS energies the relevant resonances are mainly A and N* resonances. If one goes to higher 
incident energies, e.g. to AGS energies ( about 10 GeV/nucleon ), one must take into account 
contributions of heavier mesons and resonances |34] . For the baryon-baryon channels leading 
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to K + or K mesons, we have adopted theoretical cross sections based on the one-boson- 



exchange model. These cross sections agree well with experiments ||36|| . The K + production 
cross section from the pion-induced channel has been evaluated within a resonance model 
||37|| , and a good agreement has been found between theory and data. The K~ production 
cross section from pion-baryon collisions is a parametrization of the experimental data 



In a previous study, it was shown that at SIS energies pion-baryon reactions can produce 
even more K + mesons than baryon-baryon reactions [ 38| . With an increasing number of 
participating nucleons, the K + multiplicity originating from the 7r-induced channel increases 
much faster than that from the baryon-baryon channel. Hence this channel is essential in 
order to reproduce the mass dependence of the K + production observed in experiments 
|39| . Furthermore, the pion-baryon reaction is important if one wants to understand the 



experimental observation of a non-isotropic K + angular distribution in symmetric heavy ion 



reactions at subthreshold beam energies j40|. This is due to the fact that the pion-baryon 
reaction of Eq.(3) has a substantial P-wave component, which results in a large anisotropy. 
The P-wave contribution turned out to be essential for the understanding of the observed 
K + angular distribution [H|. 



C. Kaon reabsorption and rescattering 

Once produced a K + meson can be hardly absorbed inside the medium because of 
strangeness conservation. However, a K~ meson can be easily destroyed through the reaction 
K~ + N — ► 7i + Y. The inverse process also contributes partially to the K~ production. We 
used in our calculation a parametrization of the experimental cross section for these reactions 
f36fl . The K~ absorption is rather strong. <Jk-n->itY is about 50 mb at a laboratory K~ 
momentum of 0.2 GeV/c and increases rapidly with decreasing momenta. 

Both K + and K~ mesons may scatter elastically with nucleons as propagating through 
the nuclear medium. The i^ + -nucleon interaction is relatively weak ( <Jkn~ 10 mb ) ||42|| . 
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However, in the initial stage of the nucleus-nucleus collisions at SIS energies, one can expect 
nuclear matter densites up to three times saturation density ||33|| , which considerably shortens 
the mean free path of the kaons (A= (cr-p)^ 1 ). This may lead to observable effects on the 
kaon spectra pi| , ^3| . Compared to K + mesons, K~ mesons have a much larger elastic cross 
section, which is about 50 mb at a laboratory K~ momentum of 0.2 GeV/c and becomes 
even larger with decreasing momenta. For the elastic K + and K~ cross sections we use the 



parametrization given in Ref. |36 



D. Kaon Coulomb interaction 



In Ref. [0 it was shown that the Coulomb interaction between the K + mesons and the 
nuclear medium has a visible effect on the kaon dynamics in heavy systems like Au + Au. 
Thus the Coulomb interaction of both the K + and K~ mesons is taken into account. The 
Coulomb field acting on the kaons is generated by protons, charged resonances and pions. 
We should notice that the multiplicity of negative charged pions has been observed to exceed 



greatly the value for positive charged pions in large systems such as Au + Au |f22| , |32| . |33 
This leads to a non-vanishing contribution of the pions to the charge density. Thus the total 
Coulomb interaction is given by the sum over all charged particles, i.e. 

jjCou _ Zi k Zi n e2 _|_ Z ik Z ir e 2 Z ik Z iir e 2 

■ ^ I /Vi . ivi . I ■ ^ I /vi . iy . I • ^ I ty . ty . I 

i n I ' *fc ' *n I i T I ' «fc ' tr \ in I ' tk 1 In \ 

where i n , i r and i n denote nucleons, resonances (A and N*) and pions, respectively. 

E. Kaon mean field of the strong interaction 

The strong interaction of kaons with nuclear matter is expected to play also an important 
role in modifying the propagation of K + and K~ mesons. Possible effects of the kaon-nuclear 
matter strong interaction have attracted special interest |T7|- [T9|j35|j3m . The well-established 



formalism for studying kaon properties in nuclear matter is chiral perturbation theory. Based 
on a SU(3)lxSU(3)j? chiral Lagrangian Kaplan and Nelson [f|4J] constructed an effective 



meson-baryon Lagrangian, whose coefficients had been determined from experimental mea- 
surements, however, still with an uncertainty of about 30%. Applying this Lagrangian to 
the KN interaction on the mean field level, one can obtain the following dispersion relation 



for a kaon or antikaon in nuclear matter Wj 



j — o i 2 n ^ KN i / 3ps v„, 3 p B , v 



/-2 > I I ^ KN i / , 21 3p B 



where denotes the bare kaon mass, and f w ~ 93 MeV is the pion decay constant, ps 
and ps are the baryon density and the scalar density, respectively. The kaon-nucleon sigma 
term Hkn depends on the nucleon strangeness content and the strange quark mass. A value 



of Ej^AT = 350 MeV has been used in the literature |i7|j36[| , however, recent lattice QCD 



calculations imply a higher value of Y>kn = 450 ± 30 MeV |^5|]. While the Y^kn term 



provides a scalar attraction for both K + and K~ mesons, the term proportional to ps gives 
rise to a vector potential, which is repulsive for K + mesons and attractive for K~ mesons due 
to the G-parity transformation. The expressions Eq.(6) and Eq.(7) for the in- medium kaon 
potentials include only the S-wave interaction. There is also a P-wave contribution of kaons 
with the nuclear medium by hyperon-nucleon hole excitations. This P-wave interaction 
changes signs from K + to K~ mesons and is in its absolute value about twice as large for 
K~ than for K + mesons. This P-wave interaction has till now been neglected in all previous 
works on kaons in nuclei. We will neglect it here also. But future investigations should 



study, if it is important. The strong potential for a kaon is in this approach defined as [46 



U** = u> K -y/mj c +p. (8) 

Considering the mean field of the strong interaction and the Coulomb field, the total Hamil- 
tonian for the kaons reads 

Hk = £ ( M + PL + K ou + u t) > ( 9 ) 



There exist experimental constraints on the kaon mean field. As indicated above, the 
K + N interaction is relatively weak compared to other hadron-nucleon interactions. There- 
fore, the impulse approximation should be justified to determine the in-medium K + potential 
at low densities from experimental free KN scattering data. This yields a repulsive K + po- 
tential of about 30 MeV at saturation density (p = 0.16 fm~ 3 ), if one adopts an isospin 



averaged free KN scattering length clkn = -0.255 fm ||47|| . On the other hand, experiments of 
kaonic atoms can provide direct information on the K~ mean field around po- An attractive 
i^-nucleus potential of U K - = -200 ± 20 MeV had been found in 56 Ni p|. With E KN = 



350 MeV, one obtains from Eq.6-8 a K + potential of about +7 MeV and a K~ one of -100 
MeV at po- Both values are quite different from their respective empirical findings. High 
order corrections and P-wave contributions to the mean field approximation may be respon- 
sible for this discrepancy. In Ref. [f45[,[49| it was argued that the so-called range term, which 
is of the same order in chiral perturbation theory as the Y>kn term, has also contributions. 
In addition, the presence of the nuclear medium can reduce the pion decay constant f n by 
decreasing the quark condensates. Taking the range term and the in-medium reduction of 
f n into account, one obtains at po about +28 MeV for the K + and -175 MeV for the K~ 
potential with a value of T*kn = 450 MeV. Thus the strength of the potentials is now roughly 
consistent with the respective empirical values. 

The i^~-nuclear matter interaction is complicated by the A(1405) resonance in the K~- 
nucleon channel, which is commonly considered to be an unstable KN 1 = bound state 
about 27 MeV below the K~p threshold. This state also appears as a resonance in the 7rS 
channel. The A(1405) is responsible for the repulsive isospin-averaged KN scattering length 
in free space ||48|| . Since the KN interaction is strong, the impulse approximation should 



not be good. A coupled channel description of the KN and 7r-hyperon channels is necessary. 
Such studies including the A(1405) have been performed by Weise et al. |5T],[52|. They show 



that the A (1405) state dissolves in a nuclear environment due to Pauli blocking. Therefore, 
although the A (1405) dominates the threshold KN dynamics at very low densities (p < 



S 



0.3po) ; it is of minor relevance for the i^~-nuclear matter interaction at densities p > po- 
These calculations pT| , |52f have been performed with free kaons. A selfconsistent inclusion 
of the kaon-nucleus potential (in medium effects) may change these conclusions and should 
in the future be considered. 



For the K + mesons, we use the potential suggested in Ref. |i5,IDfl, since it agrees with the 
impulse approximation. For K~ mesons, a potential of -200 MeV at po seems to be too deep 
to describe the experimental data of K~ production from heavy ion reactions PB| , |SD| . Thus 
one usually uses a K~ potential of about -100 MeV at normal nuclear matter density [5B,50|. 



We use a K~ potential derived from Eq.7-8 with T*kn = 350 MeV. The additional corrections 
from P-wave contributions, the range term and the in-medium modification of the pion decay 
constant are thereby neglected. Such a K~ potential is more or less the same as used by other 
authors |36|j50|] , and agrees reasonably well with coupled channel studies ||5Ti|52| . We will 
show later that the conclusions of the present work do not depend on the particular choice of 
the K~ potential. However, there seem to exist some discrepancies between the K~ potential 
extracted from the kaonic atoms and the one required for understanding K~ production from 
heavy ion reactions |36|j50|| . Future studies should shed light on this question. In Fig. 1 the 
corresponding potentials for K + and K~ mesons used in the present work are shown for zero 
momentum. We note that the scalar density p$ is model dependent. We followed a method 
used in the RBUU ( Relativistic Boltzmann-Uehling-Uhlenbeck ) model to calculate ps as 
a function of ps |53| . 



The kaon production thresholds might also be changed in the nuclear environment by 
the strong mean field U . ( We note that the Coulomb field JJ Cou has no net effect on 



the threshold due to charge conservation. ) The KaoS Collaboration [35] has observed an 
enhanced K~ production in nucleus-nucleus collisions as compared to free nucleon-nucleon 
collisions, which could be attributed to an attractive K~ mean field. Thus, we have taken 



into account the reduction of the K production threshold in a way similar as in Ref. [36 
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The K + mean field is much weaker than the K one. About medium modifications of 



the K + production threshold exist different opinions in the literature p6| , ^5| , |5^| . If such an 
effect exists, it should mainly affect the multiplicity of the kaons, rather than their collective 
motion. Therefore, we have not included any medium modification of the K + production 
threshold in the present work, since the K + mean field is much weaker than that for the 
K~ mesons. One can expect that the threshold is also much less affected. 

The QMD model described above has been shown to be able to reproduce very well the 
experimental measurements of the multiplicities and transverse flow of nucleons, pions and 
kaons |^,[3^,[|l],[33|] . In this work, we will apply the same model to other kaon collective 



flows (e.g. radial kaon flow). 



III. K+ RADIAL FLOW 

A. Transverse mass spectra 

One way to obtain information on the collective motion in heavy ion reactions is to 
investigate particle multiplicities as a function of the transverse mass ( m t = \J uiq + P% )■ 
In Fig. 2 we show the transverse mass spectrum of K + mesons emitted at midrapidity ( 
-0A<y cm /y pro j<0A ) in the reaction Au + Au at an incident energy of 1 GeV/nucleon and 
an impact parameter of 3 and 5 fm, respectively. A striking finding of the present work 
is that the kaon spectrum exhibits as a function of the transverse mass a "shoulder-arm" 
shape, which deviates from a pure thermal picture. This becomes evident by comparing the 
kaon spectrum with a Boltzmann fit to the high energy part of the spectrum ( also shown 
in Fig. 2 ). In a pure thermal picture the kaons are produced from an equilibrated thermal 
source. The transverse mass spectra should thus satisfy a Boltzmann distribution, which is 
more or less a straight line if plotted logarithmically. A similar deviation from the thermal 
behavior had already been found for nucleons and composite particles ( deuterons, heliums 
etc. ), which was explained by introducing a collective expansion of the fireball, i.e. the 
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nuclear matter radial flow [55|-p7|. This expanding source picture has also been successfully 
applied to pion and kaon spectra at AGS energies ( about 10 GeV/nucleon ) f58]j . This 
seems to be reasonable since at such high beam energies hadron-hadron collisions occurring 
in the late expansion phase of the reaction are still energetic enough for kaon production. 
The kaons produced in the expanding nuclear matter should reflect a collective component 
transfered from the nucleons to kaons. However, it is not clear if this mechanism still works 
at subthreshold energies. 

In Fig. 3 calculations are presented where the K + potential U str and the Coulomb 
interaction JJ Cou have been dropped. Now one sees only effects of the hadron-hadron col- 
lisions, which produce the kaons. ( The i^ + -nucleon elastic scattering is also included in 
the simplified calculations, however, its influence is fairly weak. ) The resultant kaon trans- 
verse mass spectrum with the same rapidity cut as in Fig. 2 is in exact agreement with a 
Boltzmann distribution. Thus, we can deduce that the "shoulder-arm" structure observed 
in Fig. 2 is caused by the kaon mean field and the Coulomb interaction, rather than by a 
collective expansion of the kaon sources. As pointed out in Ref . |18j , the mean field and not 



the Coulomb interaction plays the dominant role in modifying the kaon motion. The fact 
that the potential and not an expanding K + source is responsible for the "shoulder-arm" 
structure can be connected with the subthreshold beam energy. Only the hadron-hadron 
collisions during the early compression phase of the reaction have enough energy to produce 
a kaon. As the fireball starts to expand coherently in a later stage, the kaon production 
ceases to stop since the energetic collisions have been exhausted. Consequently, at a sub- 
threshold beam energy no " shoulder-arm" spectrum for the midrapidity kaons can arise from 
the velocities of the sources. On the other hand, the K + mean field is repulsive with a value 
of about 28 MeV at saturation density ( po ) an d increases with the density. The kaons will 
thus inevitably experience an acceleration due to the repulsive potential as they propagate 
outwards from the overlap region. This may give rise to a collective motion of the kaons 
in the radial direction, which can lead to a deviation from the thermal picture. Therefore, 
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the observed " shoulder- arm" structure in the kaon spectrum appears to be a signature of a 
genuine collective motion of the kaons, rather than only a result of the nuclear matter radial 
flow. 



B. Radial flow 

In order to quantitatively extract the collective component from the kaon spectrum, we 
fit the QMD results by incorporating a common radial velocity of the K + mesons to the 
standard Boltzmann distribution. The corresponding distribution reads 



d 3 N n m , „ E 2 E 2 JhE -aT) 2 -m 2 



d<pdym t dm t p 2 p 2 p 

where E = m t coshy, p = \/p 2 + m t 2 sinh 2 y, a = j(3p/T, 7 = (1 — f3 2 )~ l l 2 . T is the 
temperature of the thermal source, while the parameter (3 = v/c is a common radial velocity 
of the kaons. Fig. 2 shows that Eq.(lO) gives a good fit to the QMD results. The fit yields 
for b = 5 fm a temperature of T = 60 MeV and a common velocity of f3 — 0.1. A slightly 
higher temperature ( T = 62 MeV ) and larger common velocity ( f3 — 0.11 ) are obtained 
for more central collisions ( b = 3 fm ). Are the T and (3 derived from the fit a good measure 
of the thermal and collective motion of the kaons? As we mentioned above, the thermal 
feature of the kaon spectra described by T reflects the temperature of the kaon sources, while 
the collective component - given by (3 - is due to the acceleration of the kaons by the mean 
field. Thus the extracted temperature describes the thermal motion without the influence 
of the kaon potential ( see Fig. 3 ). On the other hand, the collective motion (3 can be 
quantitatively described by the enhancement of the kaon momentum by the repulsive K + 
potential. The simulations shown in Fig. 3 give temperatures of T = 64 MeV for b = 5 fm 
and T = 66 MeV for b = 3 fm, respectively. They are quite comparable to the temperatures 
T extracted from the fit to the full calculations (Fig. 2). The average transverse mass of the 
kaons is increasing from <m t > = 566 MeV to 594 MeV for b = 5 fm as one changes from 
the simulations without to those including the kaon potential. For b = 3 fm it increases 
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from <m t > = 571 MeV to 600 MeV. ( We note that the Coulomb interaction contributes 
less than 30% to the enhancements. ) This corresponds to a common radial velocity boost 
of the kaons by (3 = 0.102 and (3 = 0.113 for b = 5 fm and 3 fm, respectively. The two values 
are very close to the (3 parameters obtained by the fit. This indicates that it is possible to 
separate the thermal and the collective components of the kaon spectrum by the present 
analysis. 

At an incident energy of 1 GeV/nucleon, the energy available for kaon production is quite 
limited. Consequently, the kaons have also a limited thermal motion. The average radial 
velocity of the kaons due to thermal motion is about 0.488 ( in unit of the light velocity c 
) for Au + Au, while the flow velocity is about (3 = 0.1 as mentioned above. The collective 
motion is more than 20% of the thermal motion, so that it is not obscured by the latter, 
and can be observed as a clear signal. 

The radial collective motion of the K + mesons may provide a new opportunity for inves- 
tigation of the kaon mean field in nuclear matter. As stronger the kaon potential, as more 
significant is the radial kaon flow. Thus the kaon radial flow can be used as a meter to deter- 
mine the kaon in-medium potential. We note that an analysis of the transverse momentum 
of the kaons as a function of the rapidity (so-called kaon transverse flow) might give also 
information of the kaon potential. Without the kaon mean field, the transverse kaon flow 
follows the flow patterns of the baryons, and exhibits a S-shape. This has been found in 
both RBUU calculations [T7| and previous studies using the QMD model [ 18 1 . The S shape 



disappears and a zero flow can be observed if one takes into account the kaon potential as 



found in Ref. (TMll- 



IV. K~ "VIRTUAL" RADIAL FLOW 

In contrast to the K + mesons studied in the preceding section, many other hadrons 
such as K~ , A, S etc. are supposed to feel an attractive mean field in nuclear matter. In 
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addition they may be strongly absorbed or scattered in the medium. If these hadrons can 
also develop a radial collective motion in heavy ion reactions is not clear. Therefore, it is 
of high interest to investigate the possible formation of radial flow of these hadrons. As a 
first example, we consider in the present work K~ mesons which are subjected to all the 
final-state interactions mentioned above. 

In Fig. 4 the K~ transverse mass spectrum at midrapidity is shown for the Au + Au 
reaction at a single impact parameter b = 5 fm and for the Ni + Ni reaction by integrating 
over impact parameters 0< b <4 fm. The considered beam energies ( 1.8 GeV/nucleon for 
the Au + Au reaction and 1.93 GeV/nucleon for the Ni + Ni reaction ) are below the free 
NN K~ threshold. We analysed the K~ mesons from a rapidity interval ( -0.3 < y c .m./y P roj 
< 0.3 ) normalized to projectile rapidity y pro j, which looks narrower than the one for K + 
mesons ( -0.4 < y c .m./y P roj < 0.4 ). In fact, the size of the intervals in units of rapidity 
is very close to each other, since we consider for the K~ production higher beam energies 
than for the K + production due to the different thresholds ( 2.5 GeV compared to 1.58 
GeV ). The rapidity of the K + or K~ mesons is explicitly included in the equations fitting 
the kaon spectra ( see Eq.10 for the K + case and Eq.ll to appear later in this section 
for the K~ case ). Therefore, small differences between the rapidity intervals for the K + 
and K~ mesons should not influence the conclusions of this study. Also shown in Fig. 4 
is a Boltzmann fit to the high energy part of the spectrum ( mr - %- > 0.1 GeV ). 
One can clearly observe a "concave" structure in the K~ spectrum again deviating from a 
thermal Boltzmann distribution. This structure is, however, completely different from the 
"shoulder-arm" shape found for K + mesons. 

The "concave" structure is a consequence of the attractive K~ mean field. Without the 
mean field this structure vanishes. ( The present work also includes the Coulomb interaction 
of the K~ mesons. However, the Coulomb effects are found to be negligible compared to 
those originating from the strong potential. ) Fig. 5 demonstrates that the K~ spectra 
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obtained without the mean field and the Coulomb interaction ( but including rescattering 
and reabsorption ) are in good agreement with a Boltzmann distribution. At a subthreshold 
beam energy hadron-hadron collisions, i.e. sources for the K~ mesons, can not induce this 
"concave" structure since the K~'s are mainly produced before the fireball starts to expand. 
We already discussed this point in Section III. It is important to realize that also reabsorp- 
tion and rescattering do not lead to the observed " concave" structure. The transverse mass 
spectrum of the midrapidity K~ mesons is determined primarily by the interaction of the 
K~ mesons with the fireball. If one neglects all the final-state interactions, the primordial 
K~ spectra should satisfy according to Fermi's argument J28| a pure Boltzmann distribution. 
On the other side, the fireball baryons show predominantly a single-temperature structure 
before the formation of the nuclear matter radial flow. We know that both rescattering and 
reabsorption are in nature stochastic processes which tend to thermalize a narrow distribu- 
tion. One can hardly expect that a stochastic interaction between two thermal distributions 
can change the thermal distributions to a two-temperature structure. Thus one observes 
in Fig. 5 a Boltzmann distribution of the K~ mesons even though rescattering and reab- 
sorption have been included. However, we should note that rescattering and reabsorption of 
hadrons by spectators can induce non-thermal features since the spectators are collectively 
deflected sidewards in the reaction plan ( nuclear matter bounce-off ). Absorption of pions 
by the spectator baryons is responsible for a pion transverse flow which is anticorrelated to 
the nucleon flow ||14j-|T6|| . Rescattering and reabsorption by the spectators can also give rise 
to an azimuthal anisotropy in the hadron emission (Squeeze-out) [f20|-|22||. 

In Section III we have shown that the K + radial flow can be described very well by 
boosting a Boltzmann distribution with a common K + radial velocity directed away from 
the fireball as a result of the repulsive K + mean field. Here we try to understand the 
K~ spectra with a similar boost, however, now with a common K~ radial velocity tending 
towards the center of the fireball since the K~ mesons feel an attractive mean field. In 
analogy to Eq.(lO) the K~ distribution reads 
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d 3 N <je, , E 2 E 2 JhE + aT) 2 -m 2 



d4>dym t dmt p 2 p 2 p 

T is the primordial X - temperature, /? is the common K~ velocity. Fig. 4 shows that 
Eq.(ll) provides a good fit to the K~ spectrum obtained from the full calculation. The 
adjusted temperature parameters ( T = 83 MeV and 76 MeV for Au + Au and Ni + Ni, 
respectively ) are again close to these values deduced from a Boltzmann fit to the QMD 
results without the mean field ( T = 87 MeV and 82 MeV for the two reactions, respectively 
). This indicates that the fitting parameter T in Eq.(ll) is a good measure for the thermal 
motion of the K~ mesons. The average transverse mass of the K~ mesons is found to be 
reduced by the mean field from < > = 596 MeV to 575 MeV for Au + Au, and from 
585 MeV to 569 MeV for Ni + Ni, which correspond to a radial decelerating boost of f3 — 
0.083 and 0.067, respectively. Both (3 values agree with that from the fit ( f3 — 0.085 for Au 
+ Au and 0.066 for Ni + Ni ). Therefore, the non-thermal feature of the K~ mesons can be 
described well by a radial velocity (3 extracted from the fit. 

The deviation of the K~ transverse mass spectra from a pure thermal distribution can 
therefore be understood in terms of a common K~ radial collective motion, which arises 
from the deceleration by the attractive mean field. We call this type of collective motion 
"virtual" flow in order to distinguish it from the K + radial flow driven by a repulsive mean 
field. Such a virtual flow is characterized by a "concave" structure in the transverse mass 
spectrum, while a radial flow shows up in a " shoulder- arm" shape of the spectrum. 

Rescattering and reabsorption alone can not induce the "concave" feature in the K~ 
spectra. However, a stochastic process may interfere with a collective motion as mentioned 
above. Reabsorption plays for K~ mesons a larger role than for K + mesons. K~ rescattering 
on the fireball baryons may increase the energy of K~ mesons. This will cancel partially the 
virtual flow caused by the mean field. Therefore, we observe in the present work a relative 
small reduction of the K~ transverse mass ( A< itlt > = - 21 MeV and - 16 MeV for Au + 
Au and Ni + Ni, respectively. ) although we have adopted a mean field for the i^~-nuclear 
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matter interaction of about -100 MeV at normal nuclear matter density. 



A " concave" structure had been observed experimentally in the transverse mass spectra 
for K + and K~ mesons in Si + Pb at 14.6 (GeV/c)/nucleon |59| and Au + Au at 10.8 
(GeV/c)/nucleon |Kj as well as for pions in many reactions and at various beam energies 
pT| . In all the cases the beam energies are much higher than the corresponding production 
thresholds for free NN collisions. Considerable numbers of kaons or pions are produced even 
after the fireball starts to expand. Thus the fireball expansion can independently cause non- 
trivial collective features in the particle spectra ]57| , |55| . Resonance decay producing pions 
has also a significant contribution, which had been found to be able to induce a "concave" 
structure in pion transverse mass spectra p0j|3ll , |62i . It is worth noting that some explana- 



tions for these observations have also used an attractive pion or kaon in-medium mean field 
However, it seems that there exist some differences between these observations and our 
present work where the subthreshold beam energy plays an essential role. 

Finally we want to add some remarks concerning the kaon potentials used in this work. 
For the description of the K + radial flow a potential has been used which agrees with 
the free KN scattering data. Based on the same potential, we were able to reproduce the 
experimental K + transverse flow data in a previous work [|1^]. It is also unlikely for resonance 
decay to modify the " shoulder-arm" shape of the K + spectra, since, at the incident energies 
studied here, very few Phi mesons and antibaryons can be created. The K~ potential applied 
is less attractive than that extracted from kaonic atoms. Since the deceleration of the K~ 
mesons by the attractive mean field is responsible for the "concave structure" of the K~ 
spectra, a more attractive potential should even favour a more pronounced "concave" shape. 
The A(1405) might also play a role in leading to a "concave" structure of the K~ spectra 
just analogous to the role of the A resonances for pions. However, the formation of the 
A(1405) is reduced at high densities. Since the K~ radial flow arises from interactions with 
the dense fireball, one can expect that the A(1405) is only of minor importance. We also 
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note that our treatment of the kaon dynamics is obviously non-covariant, since e.g. the space 



component of the vector part of kaon-nucleus potential is neglected as usual ]B3] . However, 
this treatment is sufficient for the present studies, since the relative motion between the 
fireball and the midrapidity kaons is quite small at the subthreshold beam energies. The 
present conclusions are not significantly affected by a full covariant treatment of the kaon 
dynamics. 



V. SUMMARY 

In this work, we have studied possible collective motion of K + and K~ mesons emitted 
at midrapidity from heavy ion reactions at SIS energies (1-2 GeV/nucleon) by comparing 
with a thermal picture. The main findings are: 

(1) . A collective flow of K + mesons in the radial direction should exist in central nucleus- 
nucleus collisions, which can be identified by a characteristic "shoulder-arm" shape of the 
transverse mass spectrum of the midrapidity kaons. The K + radial flow is caused by the 
repulsive K + mean field. 

(2) . K~ mesons also exhibit a collective motion in the radial direction in heavy ion 
reactions at subthreshold beam energies even though they suffer strong reabsorption and 
rescattering in the nuclear medium. The collective motion leads to a characteristic "concave" 
structure in the transverse mass spectrum of K~ mesons emitted at midrapidity. This 
distinguishes the K~ collective motion from that of K + mesons. The K~ collective motion 
is induced by an attractive mean field. 

(3) . Since the radial flow of both K + and K~ mesons arises from their respective in- 
medium mean field, the radial flow may act as a novel tool for investigation of possible 
modifications of kaon properties in dense nuclear matter. 
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Figure Captions 

Fig. 1. The in- medium potentials of K + and K~ mesons in nuclear matter used in this 
work. 

Fig. 2. Transverse mass spectra of the K + mesons emitted at midrapidity ( - 
0A<y cm /y pro j<0A ). The histograms are results of the QMD simulations for Au+Au 
reactions at 1 GeV/nucleon and at two different impact parameters b = 3, 5 fm. The 
calculations are performed including the full in-medium kaon dynamics. The lines are fits 
to the transport calculations. In the upper panel the result of a pure thermal Boltzmann fit 
to the high energy part of the spectra ( m t -m > 0.15 GeV ) is shown. In the lower panel 
the fit according to Eq.(10) is presented where one assumes a K + radial flow in addition to 
the thermal motion. 

Fig. 3. Transverse mass spectra of the K + mesons emitted at midrapidity ( - 
0A<y c , m Jy pro j<0A ) for the same reactions as in Fig. 2. The histograms present results of 
the transport calculations without the K + mean field and Coulomb potential. Now a pure 
thermal fit (lines ) is sufficient to reproduce the spectra. 

Fig. 4. Transverse mass spectra of the K~ mesons emitted at midrapidity ( -0.3 < 
Yc.m./Yproj < 0.3 ) in a semicentral ( b = 5 fm ) Au + Au reaction at 1.8 GeV/nucleon 
and in a Ni + Ni reaction at 1.93 GeV/nucleon. The Ni + Ni reaction is averaged over the 
impact parameter region < b < 4 fm. The histograms are the results of the transport 
calculations including the full K~ in-medium dynamics and the lines are the corresponding 
fits to the spectra. In the upper panel the result of a pure thermal Boltzmann fit to the 
high energy part of the spectra ( m t -m > 0.1 GeV ) is shown. In the lower panel the fit 
according to Eq.(ll) is presented where one assumes a K~ "virtual" radial flow in addition 
to the thermal motion. 

Fig. 5. Transverse mass spectra of the K~ mesons emitted at midrapidity ( -0.3 < 
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Yc.m./Yproj < 0.3 ) for the same reactions as in Fig. 4. The histograms present results of 
the transport calculations without the K~ mean field and Coulomb potential. Now a pure 
thermal fit (lines ) is sufficient to reproduce the spectra. 
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